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Abstract — Localization of moving ground targets us¬ 
ing acoustic signals obtained by a passive sensor net¬ 
work, made up of acoustic sensor arrays on the ground, 
is a difficult problem as the signals are contaminated 
by wind noise and hampered by road conditions and 
multipath, etc., and are generally not deterministic. It 
becomes even more challenging when some of the ve¬ 
hicles are wheeled (e.g., cars) and some are tracked 
(e.g., tanks), and are closely spaced. In such cases the 
stronger acoustic signals from the tracked vehicles can 
mask those from the wheeled vehicles, leading to poor 
detection of such targets. A novel detection scheme is 
presented, according to which the direction of arrival 
(DoA) angle estimates of emitters are obtained by each 
sensor array using real data. The full position esti¬ 
mates of targets, obtained following the association of 
the DoA angle estimates of the same target from at 
least three sensor arrays, are used for target tracking. 
However, because of the particular challenges encoun¬ 
tered in multiple ground vehicle scenarios, this associ¬ 
ation is not always reliable and thus, target tracking 
using kinematic (DoA angle) measurements only is dif¬ 
ficult and it can lead to lost tracks. In this paper we 
propose a new feature-aided multidimensional assign¬ 
ment algorithm, to augment the existing assignment al¬ 
gorithms which use only kinematic measurements, to 
improve the association performance, especially in the 
case of wheeled vehicles. We present a novel frequency 
domain feature extraction technique by implementing a 
statistical characterization of the features, in order to 
enhance the accuracy of data association. The feature 
and DoA angle measurements are used simultaneously, 
via a joint likelihood function, in a multidimensional 
assignment (MDA) to localize targets. 

Keywords: Multidimensional assignment (MDA), feature- 
aided MDA, min. variance distortionless response (MVDR). 
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1 Introduction 

Ground vehicle tracking using acoustic data obtained 
from passive sensor networks is a very challenging prob¬ 
lem. Passive acoustic sensors are gaining in popularity 
because of their low cost, the ease of deployment, and 
the fact that they can be deployed on the ground. As 
passive sensors do not emit their own signals like active 
sensors, there is no danger of being detected. Passive 
acoustic sensors can be used in battlefield monitoring 
as well as in surveillance for civilian applications. 

In single target scenarios with active sensors, kine¬ 
matic measurements (such as range, bearing, etc.) can 
be obtained, which can be used to estimate the trajec¬ 
tory of targets [2]. When a network of passive sensors 
is used, however, the range of a target can be obtained 
only after associating the direction of arrival (DoA or 
line of sight — LoS) angle estimates obtained by at 
least three sensors. However, “ghosting” can occur, 
especially in the case of road convoys [1]. Data asso¬ 
ciation is difficult when the targets stay close together 
over an extended period of time (as acoustic signals 
from some targets can fade and then re-appear or can 
be masked by stronger signals from other targets), be¬ 
cause one has to associate the DoA estimates from the 
same target to obtain its position. Tracking using 
DoA angle measurements can be done using various 
classification techniques [3, 5], if the identities of signal 
sources (targets) are known. However, this is an unrea¬ 
sonable assumption in multiple ground target acoustic 
sensor network scenarios, due to the various challenges 
discussed earlier. Feature-aided tracking (FAT) tech¬ 
niques are a new and rapidly developing research area, 
as they exploit certain properties of the signals, called 
features, in addition to the kinematic (DoA) measure¬ 
ments, to alleviate the difficulties encountered in con¬ 
ventional target tracking with kinematic measurements 
only. 

Many of the algorithms in the literature assume the 
presence of a single target [5, 6] and use statistical pa¬ 
rameters, namely the mean and variance of several har- 
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monies of the fundamental engine firing rate of each 
target for identification and classification. When mul¬ 
tiple targets are present within the surveillance region 
of acoustic sensor arrays, the measurements no longer 
exhibit the same statistics as they did for the individ¬ 
ual targets. As a result, separate algorithms have to 
be developed for localizing ground vehicles in a con¬ 
voy. Feature-aided data association techniques em¬ 
ploy a feature-augmented measurement set, as com¬ 
pared to just kinematic measurements. Possible fea¬ 
tures could be the acoustic signature of targets from 
the power spectral density (PSD) of the acoustic sig¬ 
nals. They can be used to augment the kinematic 
measurements obtained by each sensor to form an aug¬ 
mented measurement set. The major challenge is the 
development of reliable models to extract and charac¬ 
terize features. Section 2 describes the generation of the 
PSD of the signals received by the passive sensors using 
the minimum variance distortionless response (MVDR) 
technique [9], and the detection of DoAs of signals us¬ 
ing the shape of the PSD. Section 3 describes a novel 
algorithm used to extract features from the PSD and 
their statistical characterization. Section 4 introduces 
and describes the target localization problem, i.e., esti¬ 
mation of full position (composite measurements) from 
DoA angles augmented by features. The composite 
measurements can be assigned to tracks using a dy¬ 
namic assignment algorithm. Section 5 describes the 
real scenario and Section 6 provides the target localiza¬ 
tion performance comparison results between the con¬ 
ventional cost and feature-augmented cost based mul¬ 
tidimensional assignment algorithms. 

2 DoA Detection via Second- 
Order Derivative Based 
Thresholding 

Circular sensor arrays made up of M microphones 
arranged equi-distantly, are employed in a passive sen¬ 
sor network on the ground to listen to a convoy of ve¬ 
hicles. Assuming that the convoy is at a sufficiently 
large distance from the sensor array, so that the re¬ 
ceived signals from the targets can be approximated by 
a planar wavefront, various wideband beamforming al¬ 
gorithms described in [9] can be used to detect signal 
sources and estimate their LoS angles with respect to 
the known sensor array position. 

At each sampling time, an FFT is performed on the 
raw acoustic data from each microphone in an array, 
and a discrete frequency range from /i to f nb with bin 
intervals of 1 Hz is chosen for processing, where nb is the 
number of bins. The frequency bins are chosen to con¬ 
form to the typical frequencies of the acoustic signals 
emitted by engines and other moving vehicle parts. As 
a result we have a wide-band processing algorithm that 
uses data from nb frequency bins. We denote the FFT 


Power Spectral Density 



Angles (degs) 

Figure 1: Power spectral density as a function of LoS 
angles and frequency bins obtained from the MVDR 
algorithm for one sensor array. 

data used for processing at each sample by X(m,f b ), 
where to G {1,..., M} represents each microphone and 
fb G {/i,..., f nb } indicates the frequency bins. The 
MVDR algorithm provides, for each sensor array, an 
estimate of the power spectral density of the acoustic 
signals 

~( djb) a(f b )'R(f b )- 1 g(f b ) 

6 d G {6»i, ...,0 D } and f b G {/i, (1) 

and 

g(fb) = {X}’ fb V (2) 

with {X} f b denoting the column of the matrix X cor¬ 
responding to fb, V denoting the steering matrix [9], 
and 

M 

R(fb) = m E X K hYX{m, f b ) (3) 

m= 1 

The DoA angles of the acoustic signals are estimated 
by detecting peaks in the spectrum, illustrated in Fig¬ 
ure 1. Peak detection can be performed by applying 
a thresholding algorithm, using the average of the es¬ 
timated power spectrum in (1) over all the frequency 
bins corresponding to an LoS direction 6, 

, n b 

P(o) = -Y J He,h) (4) 

Ub b—1 

For each sensor array s, i s (k) DoA angle estimates are 
obtained at time k . 1 

The thresholding algorithm used in the present work 
to detect DoAs from acoustic signals received by a sen¬ 
sor array is described next. Figure 1 shows an example 
of the estimated PSD at a particular scan for a sensor 
array. DoAs are detected at each scan, using a thresh¬ 
olding algorithm applied to the PSD of the acoustic 

1 The time index k is omitted in the future for simplicity. 
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signals received by each sensor array. At each angle 6 
in the steering vector (see Ch.2, [9]), the PSD as a func¬ 
tion of the angle can be obtained by averaging over all 
the frequency bins (it is denoted by P(8) in (4)), and 
is illustrated in Figure 2(a). The first derivative of the 
average estimated power, -jgP(8), for a particular angle 
Od (where d £ {1,..., n s }, and s and k are the array and 
scan indices, respectively) is shown in Figure 2(b). A 
DoA detection is made at angle dd if a positive peak is 
detected at dd in shown in Figure 2(c). In 

the example shown in the figure, the DoA angles de¬ 
tected are 85°, 57° and 102°, arranged in decreasing 
order of their corresponding amplitudes. 



(a) (b) 



(c) 


Figure 2: The thresholding used to detect DoAs from 
a particular sensor array at a certain scan: (a) average 
power spectrum as a function of angles only, (b) the first 
derivative of the function in (a), and (c) the negative 
of the second derivative of the function in (a). 

Figure 3 shows the DoA angle estimates obtained by 
the sensor array from Figure 1. In a typical ground 
vehicle tracking scenario where the targets are moving 
in a convoy on road or off-road conditions, there are 
a variety of extraneous factors which affect the acous¬ 
tic signal, such as road conditions, sound generated by 
moving parts, wind, etc. This results in numerous false 
alarms, i.e., power detections in frequency-angle bins 
not directly due to the engine or moving parts of the 
targets. There are also missing DoA angle estimates 
(missed detections) due to signal attenuation or the 
possible masking of signals from wheeled vehicles by 
those from the tracked vehicles, especially when the 
targets are closely spaced in a convoy. As a result, the 
quality of the data association is low if just the DoA an¬ 
gle estimates are used as measurements. Therefore, the 


DOA Angle Estimates 



Figure 3: DoA angle estimates for one sensor array. 


power spectrum will be exploited to generate features 
which could enhance the accuracy of data association. 
Motivated by this, a new feature extraction procedure 
is described in the next section. 

3 Feature Extraction 

In most feature-aided tracking applications the am¬ 
plitude of the PSD peaks are used as features; how¬ 
ever, due to signal attenuation and masking in the case 
of multiple vehicle road convoys, they are not reliable 
as features. The location and the spread of the peaks 
carry more useful information about the signal source, 
compared to the amplitude, because they are not as af¬ 
fected by signal attenuation. In this paper we propose 
a statistical modeling of the distribution of the power in 
the frequency bins at each detected DoA angle and use 
a Gaussian mixture model (GMM) to extract feature 
vectors instead of just scalar features. This method en¬ 
ables the use of the location and the sharpness of the 
peaks in the power spectrum as features. 

The MVDR algorithm obtains an estimate of the 
power spectrum P(9d, fb) as a function of both fre¬ 
quency bins (fb) and LoS azimuth angles (Od), as 
shown in (1). For each sensor s, the power spectrum 
along each DoA angle 0 sis (where s £ {1,...,S} and 
i s £ {l,...,n s }), P(0 sis ), is used to declare detections 
in certain directions and to extract features correspond¬ 
ing to each detected angle at scan k. 

3.1 Fitting of a Gaussian Mixture 
Model (GMM) 

The observed data d(x) (the estimated power spec¬ 
trum in a particular DoA direction x) is modeled as 

d(x) = y(x; /3) + e(x) (5) 

where y(x\(3) is the fitted parametric model, (3 is the 
parameter vector and e(x) is the fitting error. The ob¬ 
jective is to estimate the parameters of the model such 
that the error (noise) is minimized in a statistical sense. 
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amplitudes) form a matrix of dimension n x S 



(c) 


Figure 4: Gaussian mixture model fitting to extract fea¬ 
tures: (a) feature data from sensor array 1 (DoA angle 
6 1 , peak locations: [40,15,36,49]), (b) feature data from 
array 2 (DoA angle 02, peak locations: [78,40,117,15]) 
and (c) feature data from array 3 (DoA angle d 3 , peak 
locations: [40,18,37,14]). 


The nonlinear least squares (NLLS) method is used 
to estimate the parameters of a nonlinear model - a 
GMM — used to fit the data. Gaussian mixture mod¬ 
eling is useful for peak finding applications. This model 
is given by 


/ 

40 

78 

40 


15 

40 

18 


36 

117 

37 

l 

49 

15 

14 


The peak location matrix given in (7) illustrates that, 
for the three sensor arrays in Figure 4, the peak loca¬ 
tions are not necessarily matched across the S sensors, 
i.e., a peak location is not necessarily grouped with 
other peak locations in the S — 1 lists (columns), in 
its immediate neighborhood. For example, in (7), the 
first row is [40 78 40]. The location of the peak in the 
second column (at 78 Hz) is not close to the location 
of the peaks in the first and third columns (at 40Hz). 
Hence, in order to properly group peak locations across 
lists (sensors), each peak location should be matched to 
other peak locations (in the other lists) in such a way 
that each matched iS-tuple of peak locations consists of 
peak locations which are close to each other. If a peak 
is located in such a way that it cannot be matched to 
any other peaks in the remaining lists, it is matched to 
a dummy element which indicates a missed peak detec¬ 
tion. This matching of peak locations, which leads to 
“matched” feature vectors, is described next. 


3.2 Matching the Location of Peaks via 
Assignment 


y(x;P) 



The peak location vectors, estimated by the GMM 
algorithm, corresponding to each of the detected DoA 
(6) angles 0 Hl ,0 2 i2 ,...,0 S i s are 


where f3 = [(a 1 , fi 1 , a 1 ),..., (a",/j n , a n )]' is the parame¬ 
ter vector of the GMM such that a 1 is the weight (am¬ 
plitude), n l is the location and a 1 is the width of the 
peak of component l , and n is the number of compo¬ 
nents of the GMM. The NLLS method can be used to 
estimate the parameters of the GMM which best fit the 
observed data in (5). Figure 4 illustrates the GMM fit¬ 
ting results when applied to power spectrum data for 
certain DoA detections obtained by three sensor arrays 
at a particular time. A GMM with n = 4 components 
was used to extract the location, width and the ampli¬ 
tude of the peaks from the power spectrum data. A 
larger n was found to lead to excessive uncertainty in 
the feature model. 

From the example illustrated in Figure 4, one can see 
that there are n = 4 frequency peaks from each of the 
S = 3 sensors, corresponding to three detected DoA 
angles 0\, 02 and 0 3 (subscripted by the sensor index). 
These peak locations (arranged in decreasing order of 


^1*1,... Sis ~ ^ 2 i 2 > •••’ 


/ 0 

0 

0 

\ 


Mli, 

/4 2 

• ■ ■ 0-Sis 



Mli, 

i4 2 

0Si s 


(8) 

V Plii 

/4 2 

..n 

• ■ ■ 0Sis 




where fr si is the location of the f-th component (peak) 
of the GMM (6) used to fit P(0 S i 3 ) and 0 indicates 
the dummy element 2 that signifies missed detection of 
a peak. The matching of peak locations across lists 
is done based on the solution of an assignment prob¬ 
lem formulated using costs attached to all the possible 
candidate ^-tuples of peak locations in (8). Costs are 
derived based on a likelihood ratio for each candidate 
S'-tuple. 


2 Indexed by zero, and shown as the top row. 
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3.2.1 Likelihood Function Derivations for Cost 
Calculation 

The likelihood a candidate 5-tuple of peak locations 
Mn,...,is from (8), is from the same target is 


A («(£)) = p (Mil,...,i s I «(*)) = I| 1 ~ P d{ 


1 [!-«?] 


i si 


P D{P (/4i. I «(*)) ( 9 ) 


where 




0 if ML, = 0 in (8) 

1 otherwise 


( 10 ) 


is the peak location indicator, P d j is the probability 
of detection of the frequency peak and re(t) represents 
the (unknown) true peak location, if the acoustic sig¬ 
nal is from a target t. As the identity of target t, is un¬ 
known, we use a generalized likelihood function A (k(t)) 
instead of the likelihood function given in (9), where 
the peak locations are assumed to be uncorrelated and 
distributed as follows 


m(mLJkW) =-A/'(Mii,;M*(*),(oiiJ 2 ) (if) 


vectors corresponding to the S-tuple of detected 
DoA angle estimates 61 ^, 62 ^, ■■■■> 6 si s 


[<t >lii > ^2i 2 ) 

•••) fisis} 




/ in 

r^lil 

hi 

P2i 2 

JSl 

■ hSi s 

\ 


,,0 12 
ri ii 

J 22 

lhi 2 

JS2 

■ PSi s 


(15) 

. . jlrr m 

\ n%i 

32n r n 

P'2 i 2 

is nm 

' Psi s 

) 



where pi 1 /* £ {0,mL: s } an< l 9 = 1) Each element 

Mgj in the s-th list of (8) appears exactly once in the 
s-th column of (15), while 0 appears n m — n times in 
each column of (15). It has to be noted that n m could 
vary with each 5-tuple of DoA angle estimates being 
considered, while n remains the same as it is a GMM 
fitting design parameter and is fixed. 

The illustrative example presented next corresponds 
to the same triplet (5 = 3) of DoA angle estimates 
61 , 62 , 9 3 , which gives rise to the peaks illustrated in Fig¬ 
ure 4. The triplet of peak location vectors of length 
n = 4 are as shown in (7). The triplet of feature vec¬ 
tors, each of length n m = 7, obtained after performing 
the matching of peak locations as described above is 


with 

V s ll 1 S f 

M^)= (12) 

L s =l 

and k(t) = [p\t ),..., p n {t)]', the estimated target peak 
location vector. The standard deviation of the peak 
location p si is a si , from (6). The indexing of p and 
a is augmented to indicate the peak index l. the sensor 
s it originated from, and the index i s of the detection 
from sensor s. 

The likelihood that the 5-tuple of peak locations are 
all false alarms is 


A (n(t) = 0) = p (Mi 


*(f)=o)=n 


8=1 


1 si 


VI 


(13) 

where V/ is the volume (in the frequency space) of 
sensor array s. 


3.2.2 Cost of Matching of the Peak Locations 
Across 5 lists 

The cost of matching a candidate 5-tuple of peak 
locations mm,...,* s is the negative log-likelihood ratio 
(NLLR) obtained from the likelihood functions (9) and 
(13) 


s =1 

-<^[ln (P D f V /) + ln P (/4i» I «(*)) ( 14 ) 


A multidimensional assignment algorithm is solved (see 
[1, 4, 7]), to obtain the following 5-tuple of feature 


^matched 


/ 40 

40 

40 \ 

15 

15 

14 

36 

0 

37 

49 

0 

0 

0 

78 

0 

0 

117 

0 

V 0 

0 

18 / 


4 Target Localization 
Feature-Aided MDA 


(16) 


by 


One of the most important issues in multisensor- 
multitarget (MSMT) tracking is data association. Re¬ 
cently, a class of algorithms called multidimensional as¬ 
signment (MDA) algorithms (also known as 5-D as¬ 
signment) have been developed to solve the data asso¬ 
ciation problem using an assignment approach [1, 7]. 
The present paper uses the AIDA approach to solving 
the data association problem. This approach, desig¬ 
nated as Multisensor Information Configuration III in 
[1], requires a static association of measurements to 
measurements, across sensor arrays at each time k, re¬ 
sulting in complete position estimates called “compos¬ 
ite measurements”. If the DoA angle measurements 
are augmented by features, the result will be feature- 
augmented composite measurements. 

In the feature-aided MDA problem, at time k, we 
are given 5 lists of feature-augmented measurement 
vectors from 5 passive sensors gathering data in a 
surveillance region. The position of each sensor array 
s = 1 ,...,S is denoted by x s (the sensor array posi¬ 
tions are assumed to be fixed throughout the track- 


74 



ing scenario 3 ). An S'-tuple of DoA angle measurements 
can be augmented with an S'-tuple of feature 
vectors = {<^> 1 ^,^Sj s }, where each angle 

measurement 0 S i s has its own corresponding feature 
vector (f> S i a , to form an S-tuple of feature-augmented 
measurement vectors 


where p l si (l = represents either a feature 

(detected matched peak location) or a dummy (missed 
detection of the peak in list s), as described in Section 
3.2.2. 4 

Assuming independence between measurement lists, 
(22) can be simplified to 


~ l Z h,—is> € {0,1,n s } (17) 

where i s = 0 represents the dummy measurement or 
a missed detection. Generalized likelihood ratio based 
costs are generated for assigning each such S'-tuple. 

The generalized likelihood that the S-tuple ig 
consists of feature-augmented measurement vectors 
from a target t is 

A a(Y =P( Z i 1 ,...i s \ t ) =P([ Z ii,-is,$iu...is]\(Zt,‘Pt)) 

(18) 

where x t is the ML estimate of the position of target 
t and fit is the estimated feature vector of target t. 
Assuming that the angle measurement errors and the 
feature vectors are independently distributed, we have 

a aW =p{Zi 1 ,...i s |x t )p($ n ,... is | (ft) = A%(t)A$(t) 

(19) 

where A %{t) and A®(t) represent the kinematic and 
feature generalized likelihood functions, respectively. 

The generalized likelihood that an S-tuple of kine¬ 
matic (DoA angle) measurements are from a 

target t with ML position estimate x t , is given by 

s 

A%{t)=p{Z ilt ..., ig \t) = n {[1 -Pd.] 1 -*™ 

S= 1 

•[P^p(d sia |x t )]“ fe) }(20) 


A^(t) = JJ {?>(<&,*, \<pt)Y 

s =1 


(24) 


Substituting from (23) in (24), and assuming uncorre- 
latedness between the components of the feature vector 
<t>si a , we have 

s fn m “(*«) 

A f(t) =11 l P D' e P(£i e 

a=lU=l J 

• (25) 

where Pjyi is the (nonunity) probability of detection of 
the features in list s, while 5 SI is the detection indica¬ 
tor function for the feature fjf ai . The feature vector 
component is assumed to be distributed as follows 


P(rf,i a I <Pt) =Af (p l si s -,P l {t),(a l si J 2 ) 
where 


lu\ — ^a=l Psis^si 


Z-^S= 1 °Sl 


i s G {1,..., n s } 
(26) 

(27) 


The standard deviation a l si of the peak location p l m 
is obtained from the GMM fitting in (6). 

The generalized likelihood that the S-tuple ls 
consists of feature-augmented measurement vectors 
which are all spurious or are unrelated to a target t 
is 


where u(i s ) is the DoA detection indicator function and 
Pd„ is the probability of detection of a bearing measure¬ 
ment, by sensor array s. Each DoA angle measurement 
is assumed to be normally distributed with the follow¬ 
ing pdf 


p(0 sis \x t ) = A f{9 sia -,h{x t: x s ),a 2 sie ) (21) 


where h is the measurement function. 

The generalized likelihood that an S-tuple of feature- 
vector measurements corresponding to an S- 

tuple of angle measurements are from a target 

t , with feature estimate ipt, is given by 

A%(t) = p (<£iu, (j>2i 2 , <j>Si s I Pt) (22) 

where (f> S i 3 is the feature vector (15) corresponding to 
the DoA angle measurement 9 S i s from the s-th list 

tysi, = [Psis 1 Psis > •••> Psis ] (23) 

3 The approach can be easily generalized to moving sensors 
(see, e.g., [2]). 


A^(t = 0) 


p(zr lt ...,is\t = °) 


s 

UiP I f = °)p (<£«. I pY) = 0)}“ (ls) 

S=1 


5 

n & i *=°) 

S=1 

71 m 

. 1=1 



(28) 


Assuming that the angle and feature vector element 
false alarm measurements are uniformly distributed, we 
have 


s 


AS(*=o)=n 

S=1 


1 

Vs 


71m 

n 


U=1 




u(i B ) 


(29) 


4 The length of the feature vector n m varies for each candidate 
5-tuple of feature-augmented measurement vectors, as explained 
in Section 3.2.2, additional notation is omitted for simplicity. 
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Figure 5: A passive sensor network of 4 acoustic sensor 
arrays listening to a convoy of 4 wheeled and tracked 
vehicles. 


where V s and V/ are the surveillance volumes in angle 
and frequency, respectively, of sensor array s. 

The cost of assigning the 5-tuple - s is given 

by the NLLR obtained from the generalized likelihood 
functions (19) and (29), and substitution from (20) and 
(25) 


A G A (t) Ag(t)A °(t) 

Ag(t = 0) n Ag(i = 0) 


(30) 


This cost function can be simplified to the following 
form: 


s 

8 =1 

+u(i s ) In [PD a V s p(0 sia | x t )] 

+u(u|^(l-^)ln(l-P Di ) 

U=1 

+ ^ln [P D iV/p(fx l sia \<p(t))]} (31) 

The most likely set of 5-tuples such that each feature- 
augmented measurement vector in a list is assigned 
to either other measurement vectors, or declared false, 
with the constraint that each assigned 5-tuple receives 
at most only one measurement vector from each list, is 
obtained by solving a global 5-D optimization problem 
using the cost function (31). A Lagrangian relaxation 
algorithm is used to solve this problem [1, 4, 7]. 

5 Scenario 

The algorithm developed above was exercised on real 
data obtained from a field experiment [5]. Four acoustic 
arrays were placed within a path that was traveled by 
multiple targets (see Figure 5). All the acoustic sensor 
arrays used in the experiment were circular arrays of 
microphones. 

Each sensor array has a field of view of 360 ° and 
is made up of 7 equi-distant microphones. DoA angle 
measurements are estimated at each sample time (Is), 


from the PSD of the acoustic signals received by each 
sensor array, as described in Section 2, with a mea¬ 
surement error standard deviation of 2°. The PSD 
of the acoustic signals is estimated using an MVDR 
beamformer, where a steering vector of length 360 (all 
possible angles/directions) is considered. A low-pass 
spectrum of 1-120 Hz, divided into bins of 1 Hz each, is 
used by the MVDR algorithm. 

Figure 6 shows the number of DoAs detected at each 
scan for the four sensor arrays used in the scenario. The 
measurement lists in both conventional and feature- 
aided assignment consist of the same number of DoA 
detections per scan. They only differ in whether fea¬ 
ture vectors augment (as described in Section 3) to the 
DoA angle measurements or not. One can observe from 
the figure that there are false alarms as well as missed 
detection in several scans, as there are four targets in 
this scenario and the number of detections varies from 
2 to 7. 



Figure 6: Number of DoA detections at each scan, for 
both conventional and feature-aided assignment: (a) 
from sensor array 1, (b) from array 2, (c) from array 3 
and (d) from array 4. 


For target localization using the conventional cost, 
an MDA algorithm is solved [1], with DoA measure¬ 
ment probability of detection Po a = 0.9 and sensor 
array surveillance region volume V s = 360 °. For target 
localization using the feature-augmented cost (31), an 
MDA algorithm is solved at each scan. The probability 
of detection of the features is assumed to be P D i =0.7 
and the surveillance volume in frequency is assumed to 
be V/ = 120 Hz (the width of the frequency spectrum). 

6 Results 

The real data based results shown in Figure 7 show 
the comparison between the performance of the con- 
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ventional cost and feature-augmented cost based target 
localization algorithms. The figure presents the com¬ 
plementary cumulative density function of the number 
of detected targets. The localization results obtained 
for 50 scans show that the feature-aided localization al¬ 
gorithm outperforms significantly the conventional cost 
based algorithm. The figure shows that the conven¬ 
tional cost based assignment algorithm performs poorly 
in detecting the full set of targets. This is because the 
stronger acoustic signals emitted by the tracked vehi¬ 
cles can mask those from the wheeled vehicles, lead¬ 
ing to missed detection of some DoA angle measure¬ 
ments by one or more sensor arrays. Hence, using only 
kinematic measurements for static assignment leads to 
the missed detection of some targets. 5 Augmenting the 
kinematic measurements with feature measurements al¬ 
leviates this problem by adding more information to 
the cost function, and hence, resulting in more feasi¬ 
ble solutions to the global assignment algorithm. The 
feature-aided target localization algorithm in this par¬ 
ticular example is able to localize all 4 targets, in 32 
scans out 50, while the conventional algorithm man¬ 
ages to do the same in 7 out of 50 scans. Overall, the 
estimated detection probability for a target composite 
measurement (in the absence of ground truth, calcu¬ 
lated as the total number of composite measurements 
from 50 scans knowing that there were 4 targets) is 
= 0.92 and ^ = 0.68 for the feature-aided and 
conventional schemes, respectively. 



Figure 7: Comparison of the estimated probability of 
localizing at least n targets out of 4 for the feature-aided 
assignment algorithm and the conventional assignment 
algorithm (without features). 


7 Conclusions 

In this paper a data association algorithm to localize 
multiple ground targets using acoustic signal data ob¬ 
tained by a passive sensor network, has been solved us¬ 
ing a novel feature-aided multidimensional assignment 

5 Assigning DoA angle measurements from at least 3 lists is 
needed for a detection. 


framework. A novel detection scheme has been pre¬ 
sented to detect DoAs and a new technique has been 
presented to extract feature vectors, instead of features, 
which augment their corresponding DoA angle mea¬ 
surements. An MDA algorithm is solved at each scan, 
using feature-augmented likelihood ratio based cost 
functions, to obtain composite measurements which are 
the full position estimates of targets. It is shown using 
real data obtained by a passive sensor network, that the 
feature-aided assignment algorithm significantly out¬ 
performs its conventional kinematic counterpart in the 
localization of multiple targets, where some targets are 
tracked vehicles and the others are wheeled vehicles. 
Ongoing work consists of assigning the full position es¬ 
timates of targets obtained by the feature-aided static 
assignment algorithm to tracks via a dynamic assign¬ 
ment algorithm. 
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